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Abstract: We study the correlators of a recently discovered family of BPS Wilson loops in = 4 
supersymmetric U{N) Yang-Mills theory. When the contours lie on a two-sphere in the space-time, 
we propose a closed expression that is valid for all values of the coupling constant g and for any 
rank N, by exploiting the suspected relation with two-dimensional gauge theories. We check this 
formula perturbatively at order 0{g^) for two latitude Wilson loops and we show that, in the limit 
where one of the loops shrinks to a point, logarithmic corrections in the shrinking radius are absent 
at 0{g^). This last result strongly supports the validity of our general expression and suggests the 
existence of a peculiar protected local operator arising in the OPE of the Wilson loop. At strong 
coupling we compare our result to the string dual of the AA = 4 SYM correlator in the limit of large 
separation, presenting some preliminary evidence for the agreement. 
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1. Introduction 

The supersymmetric Maldacena- Wilson ^ loops in = 4 supersymmetric Yang-Mills theory 
(SYM) were recently generalized to include a class of contours contained in an S^, which also include 
a path-dependent coupling to the scalar fields of the theory ^. A subset of those Wilson loops 
are contained in a great and their discoverers pointed out an exact solvability and a potential 
connection to QCD2 ^. These loops are given by (we consider our S"^ in hyperplane = 0) 

W = ^TtV exp j) dr [i±'Ai + €ijkX^x''M}<l>j^ (1.1) 

where x*(r) (where i = 1, . . . , 3, / = 1, . . . , 6) is a closed path on S"^, and Mj is a 3 x 6 matrix 
satisfying MM'^ = 1 and which we will take to be = 1/R (no summation implied and R is the 
S'^ radius) and all other entries zero. At the level of the vacuum expectation value (VEV) there is 
considerable evidence that^ 

where ^1 is the area on the sphere enclosed by the Wilson loop, while ^ = ^1 -|- ^2 is the total 
sphere area. To begin with, the 1/2 BPS circle (given by an equator) has been proved to be given 
by ( |1.2D [^, 1^, |8| and there are strong arguments in favour of the 1/4 BPS circle of Q (given by a 
latitude) also being captured by i^). At 0{g^), (0) was proven for general contours in p, H]- 
This result was further confirmed at 0{g^) in []lO| , |lT| . The significance of the result is that it agrees 

^L™ is the Laguerre polynomial L^(x) = 1/n! exp[a::]x~'"(d/da;)"(exp[-a;]a;""'''"). 
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with the calculation of the VEV of the Wilson loop in QCD2 on an in the zero instanton sector 
with the couplings related by^ 

9ld = -^- (1-3) 

The idea that a class of = 4 SYM Wilson loops might be exactly solvable and equivalent to 
Wilson loops in a lower dimensional theory is very attractive, and hints at a relationship between 
two very different quantum field theories. More specifically one could infer that the localization 
procedure presented in ^ could also apply to this more general class, pointing towards the existence 
of a sector of non-local topological observables in = 4 SYM. Standard field theoretical arguments 
should then suggest the presence of protected local operators arising in the OPE of the Wilson loop 
(see 1 13] for related research in this direction). 



To substantiate these ideas we need to go beyond the level of the one-point function of Wilson 



loops and consider correlators of loops. A first step in this direction was undertaken in [11 1, where a 
perturbative computation of the correlator of two latitudes at order 0{g^) was undertaken. Lacking 
a zero-instanton QCD2 result to compare to, in |jll| the generalization to S"^ of the Wu-Mandelstam- 
Leibbrandt (WML) [14, 15, prescription for QCD2 in the plane proposed in [|3|, Q was used. 



Indeed, this prescription has been recently shown to be equivalent to the zero-instanton QCD2 
result [|6|3_ 



In the present paper we derive a general formula for correlators of BPS Wilson loops with 
arbitrary contours on 5^ in terms of the multi-matrix model governing the zero instanton expansion 
of QCD2. The result is valid for any coupling constant g and for any value of N: we compute 
explicitly the matrix integral for the correlator of two loops. Our general expression survives a 
series of non-trivial tests. First of all we calculate in AA = 4 perturbation theory the correlator of 
two latitude Wilson loops at 0{g'^), finding perfect agreement with the matrix model result. Next we 



provide compact formulas for the perturbative 0{g^) contribution, generalizing the results of [1C|, 
from which a numerical evaluation can be easily performed (we will report on this point in the future 
[|l^ ). Here we prefer instead to investigate analytically the limit where one of the two latitudes 
shrinks to zero size: because our nonperturbative formula is an order by order polynomial in the 
shrinking radius, the absence of logarithmic terms is a crucial test of the matrix representation. 
We find indeed the absence of leading logarithms in the shrinking radius, a quite non-trivial result, 
differing dramatically from the analogous computation of non-BPS correlators [|T^ where logs are 
present. 

Interestingly, by analyzing the OPE of the shrinking Wilson loop one can relate the absence of 
the logarithmic terms to the protection of a local operator which may be expressed as the trace of 



the square of a twisted field strength. Work [|13|| concerning super-protected local operators could 



be extended to also include this novel operator, which is based on very similar symmetries. We 
discuss this issue in section ^. 

Armed with our general result we can therefore take the large N and strong coupling limit and 
try to compare it to the AA = 4 correlator from the string side. In the limit that the two latitudes 
shrink to opposite poles on the sphere, this calculation reduces to the semi-classical exchange of 
supergravity (SUGRA) modes between the two string worldsheets describing the Wilson loops at 
strong coupling. We find that at leading order in the large-separation limit, the matrix model result 
seems to capture the exchange of the SUGRA modes dual to a certain chiral primary operator. 



^We use different conventions for the Yang-Mills actions in two and four dimensions that differ by a factor two, 
in keeping with the original references on the subject. 



disagreement was erroneously present in [pl|. 
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Other modes, dual to other protected operators present in the weak couphng OPE, should also be 
carefully included to test definitively this result at strong coupling. We find moreover an intriguing 
pattern of matching between the QCD2 result and the exchange of heavier modes dual to chiral 
primary operators of higher dimension, which seems to extend to arbitrary order in the large- 
separation expansion. We have not yet understood the meaning of this highly non-trivial pattern 
of matching. 

In this paper we present a survey of our investigations, deferring a complete analysis with all 
the relevant technical details to a future publication. 

Note added: as this manuscript was being completed |Q appeared, presenting a partial overlap 
with the results of this paper. 



2. Symmetries of the loops and of their correlators 

We start by considering = 4 SYM Wilson loops that are a special case of the general construction 
presented in ||3|, They are 1/4 BPS supersymmetric loops with the contour defined on a latitude 
of 5^, first put forward in P]. Writing the Wilson loop as 



W = -^TrPexp j, dr {ix^A^, + \x\ 6^$/) 



(2.1) 



the latitudes are given by the following closed paths on an 5^ C and on another 5^ C which 
gives the coupling to the scalar fields = 1, . . . , 4, / = 1, . . . , 6), 

= R (sin 6*0 cos r, sin 6*0 sinr, cos ^Oi 0), O^ = (— cos 60 cos r, — cos 6*0 sin r, sin 6*0, 0, 0, 0). 




Two such Wilson loops are pictured in figure |T|. The supersymme- 
tries preserved by these operators are fully described in Q, see section 
2.3.1: here we just repeat some details of that analysis which are rel- 
evant to our work. 

Under general superconformal transformations we have for the 
J\f = 4 SYM bosons 



^'jf, e, 6,<l>i = e, e = eo + x^Tm^i- 



(2.2) 



Figure 1: Two Wilson 

loops given by latitudes at Demanding that 5^W = one finds two relations 
polar angles 9^ and 9q. 

712 ei = -ri2 ei. 



i7l2 + cos 6*0 73r2(723 + r23) 



(2.3) 



It is clear that each of them reduce the supersymmetry by half, and therefore a single latitude is 1/4 
BPS. We will be mainly interested in the correlator of two such Wilson loops, as shown in figure ||. 



The first relation in (2.3) is shared between two such latitudes, whereas the second is clearly not. 
Thus two latitudes are collectively 1/8 BPS, each sharing half of their individual supersymmetry. 
The same reasoning applies of course to a collection of n latitudes, resulting always in a 1/8 BPS 
system. 



_ ^ _ 



2.1 Operator product expansion 

In the next section we will present results of a perturbative calculation of the correlator of two 
latitudes and, in particular, we will consider the limit where one of the latitudes shrinks to a point 
at the pole of the sphere. The emerging structure can be usefully understood in terms of the OPE 
and its physical meaning is quite transparent. 

The crucial observation is that, viewed from a comparably large distance, the unshrunken 
Wilson loop sees the shrunken loop as a collection of local operators [21|: the quantum behavior 



is encoded into Wilson coefficients and anomalous dimensions. The story was worked out in detail 
for two circular Wilson-Maldacena loops in Here, for the 1/4 BPS latitude, we will find that 
the relevant OPE is quite different, giving rise to novel operators which appear to have protected 
dimensions. 

When analysing the OPE, we can in fact consider the general situation of loops with arbitrary 
contours on that are generically 1/8 BPS. As noticed in Q the Wilson loop (|1.1| ) can be written 
in terms of a new gauge connection 

Ai = Ai + icijkX^ —. (2.4) 

The OPE expansion will appear particularly simple using this generalized connection^. The first 
step is to determine the classical expansion of our Wilson loops in terms of local gauge-invariant 
operators when the circuit is small. To achieve this goal we shall assume that the circuit can be 
written as follows 

x\t)=4+rx\t), (2.5) 

xq being the point about which the loop is shrinking and r a parameter that will control the limit. 
We expand the contour integral by exploiting the Fock-Schwinger gauge (x — xgYAi^x) = 0, where 
the following formula holds in terms of the new gauge curvature J^ji 

Ai{x) = / dXX{x - xoyj^ji{xo + Xix - Xq)). (2.6) 
Jo 

The leading order result is given by 

dtA{x)x' ='^Ti,{xo) i dtx\t)i{t) + 0(r3) = i*Jfe^.^.(xo)nfc(xo) + 0{r^), (2.7) 
c ^ Jc ^ 

ni{xo) being a normal vector to 5^ at the point xq, depending on xq and the contour. The ex- 
pansion could of course be extended to any given order in r, producing a series of local operators 
Oq{x) determined by the particular shape of the Wilson loop, the generalized connection Ai itself 
depending on the contour. Because these operators should share the BPS properties of the associ- 



ated Wilson loop, we obtain a practical realization of the proposal of ||l^: in particular we could 
expect that their correlation functions, when restricted to the relevant S^, are somehow protected 
from quantum corrections. This would imply severe constraints on Wilson loop correlators. Let us 
exemplify the consequences for latitude correlators (we will consider here for simplicity the SU{N) 
case) . 

In our specific example we take as our shrinking point the north pole, xq = i?(0,0, 1), while 

00 
2 



r = sin^o and x*(t) = -R(cos t, sint, tan ^). Due to the trace in the path-exponential the first 



We thank Nadav Drukker for suggesting this course of investigation to us. 
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non- vanishing contribution to the OPE is quadratic in the fields, and we get expUcitly at leading 
order 



vr^ sin 6'h 



where 

O^(xo) = Tr hi? $3 - iR^ Fi2 -R^{di^i+ 82^2) ^ 



We note a peculiar feature that makes this OPE quite different from the usual circular Wilson- 
Maldacena case [|l^: operators of classical dimension 2, 3, and 4 all couple with the same power 
of the parameter which sets the size of the shrinking latitude: the polar angle 9 (in the standard 
case the power is the classical dimension itself). Indeed the overall scale R of the S"^ is just a place 
keeper. The conformality of = 4 SYM prevents it from playing any role, and it drops out of the 
calculation of any observable. 

We notice that we can easily obtain the leading term of the two latitude correlator at order 
g'^ from the OPE ( |2.g| ), once we restore the canonical normalization for the fields. We just need to 
compute the correlation function 

(O^(xo) jdt [x\Mxi) - ie,,kx\x\ ^\xi))) = ^^^^ + \ (2.10) 

that enters in the Wick contraction. Taking the relevant color traces we get 

; — 1 = — 2sm — cos — = cos — . (2.11) 

{Wq){Wi) 8 V 2 2 y 32 2 ^ ^ 

The above result will be confirmed in the next section by the finite size correlator. Actually we can 
learn something more: the general expectation for the structure of the OPE of a shrinking Wilson 
loop is given by |jl9|, ||, |24 1 

^ l + Y.ij{9^)L'^-'Oj{x) (2.12) 



where L is the size of the shrinking loop, and Oj{x) is an operator of classical dimension J and 



quantum dimension Aj = J + Aj^ + . . .. The Wilson coefficients Cj{g'^) depend on the coupling 



constant g^ . The curious structure of the latitude OPE is a refiection of the fact that the coefficients 
0(5^) which describe the coupling of the Wilson loop to a specific operator Oj{x) are themselves 
functions of 6 [27], and can be expanded as ^j{g^,0) = Ylik^^Pid^) ^he limit 6* — > 0. This 
provides us with the general structure for the OPE of Wq 

= 1 + E ^^^9\e,) Oj{xo) = 1 + Y: ^'j\9') 0^'^' Oj{xo) = 
^•^0^ J J,k (2.13) 

= 1 + e^'+' O2{xo) + ^0^^+' 03(^0) + ^o""' O4(xo) + . . . , 

where we have dropped the scale R (to restore it replace Oj{xq) — > R'^-' Oj{xo)), and have noted 
the vanishing of ^2*^'^^ and ^3'^^ from the explicit expression of ( |2.9D . The explicit form of 02,3,4(2^0) is 
simply obtained from Oj^{xo). Actually there are multiple operators of the same classical dimension, 
so there is an extra suppressed index on the (,j{g'^,6o), Aj, and Oj(xq), which is implicitly summed 
over in ( 2.13| ). In the last line we are referring only to the operators appearing in ( |2.9| ) as these 
are the only ones present at leading order in 9q. We derive the following general relation in the 
shrinking limit 

=1 + ^o""^^' {Wi 02{xo)) + e^^^-'' {W, O3(xo)) + ^o""* {Wi 0,{xo)) + .... 

(2.14) 
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We notice that when expanded at small coupling the 9^'^ terms generically produce logarithms 
Oq'^ = Oq + A J "* 9q log ^0 + • • • if quantum corrections modify the classical dimensions. The 
quantities ^3^\ and ^^^'^ may easily be read-off in our case from (|2.9|) . Since the operators 
appearing in the explicit expression are quadratic in the fields, one has that ^2^'*, and ^4^^^ lead 
as g'^. We therefore generally expect terms of the form 17^ log to show up in the perturbative 
expansion of the correlator at order g^, in the shrinking limit. 

The presence of logarithmic corrections would be a signal that anomalous dimensions are play- 
ing a part, suggesting that the full interacting theory should be taken into account and localization 
techniques would not be sufficient in the exact computation. It would also rule out the relation 
with two-dimensional Yang-Mills that produces just polynomial dependence on at any order of 
perturbation theory, as we will see in section ^. In section |3| we show that, surprisingly, no such 
logarithmic terms appear at order g^ , supporting the matrix model proposal. This indicates that 
the composite operator 0{x)^ arising from the OPE of the BPS loops (1.1), should be protected 
- at least at the first non-trivial quantum order. In other words logarithmic divergences should 
be absent in the two-point function {0{xi) 0{x2)) ■, when X12 belong to the relevant 5^, in the 
same way as the operators defined in [|l3|. It is not difficult to show in fact that 0{x) inherits the 
BPS properties of the latitude loop, and a certain amount of supersymmetry is preserved by its 
correlators. 



3. Perturbative results on Wilson loop correlators 

t:::^^ III this section we perform a perturbative analysis up to order^ g^ for the 

connected correlator W(Ci,C2) = W{Ci,C2) - W{Ci)W{C2) of two latitudes in 
the case that the gauge group is U{N). To begin with, we shall consider the g^ 
diagram depicted in fig. |2|. [Notice that this contribution would be absent in a 
5C/(iV)theory.] 

Figure 2: Xn order to carry out the computation, we parameterize the two circuits using 

g —diagram polar coordinates 

Ci =i?(sin 6*1 cos r, sin 9i sin r, cos ^i) 

(3.1) 

C2 =i?(sin 62 cos cr, sin 92 sin a, cos 6*2) , 

and define the effective propagator AcjC'2(t, cr) connecting the two loops 

2 sin 6*1 sin 6*2 (cos (r - fj) (cos 6*1 cos 6*2 - 1) + sin 6*1 sin 6*2) 

Ac.c.(r,a) =-(Tr(^)(r)TV(^)(a))o = gvrMcos 0^ cos ^2 + cos (r - a) sin^, sin^2 - 1) ' 

(3.2) 



where A denotes the effective field iA^{x)x^^ +Qi(^\x)\x\. Then the contribution is given by 

MCi,C2)\g2 = £^ r drda Ac,c,{r,a) = ^^ (X^g^N), (3.3) 

where A is the total area of the sphere, and Ai and A2 are the areas enclosed by the two Wilson- 
loops given by 

Ai 27r(l-cos0i) .2^1 A2 27r(l + cos02) 2^2 

— = = sm — — = = cos — . (3.4) 

A Att 2 a Att 2 ^ ^ 
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Figure 3: diagrams 



At order g'^, we have to consider the diagrams in fig. ^. First, we shall consider the contribution 
Sg2_g2 due to diagram (61). Its evaluation reduces to the following integral over the circuits 

^4 /■27r 



s. 



g^-g2 / dTidT2daida2 



9 



f-2n 

/ dTidaiAc^c2in,cri] 
Jo 



a 



2 sm — cos — 
2 2 



2iV2 ^4 ■ 



(3.5) 



Next we shall consider the contribution Sg_g3 due to the two diagrams (62). The sum of the two 
diagrams yields 

Sg_g3 =^ dn j) d(TldcT2d0-3(AciC2(n,cri)Ac2C2(f72,CT3) + AciC2(n,cr2)Ac2C2(o-l,cr3) + 



+AciC2(n,CT3)Ac2C2(o-i:Cr2)) + (Ci ^ C2) 
(2sm ycos y I (sm + sm 62)-- 



A2 



2iVM4 



AiA2(AiA3 + A2^3 + 2^i^2), (3.6) 



where Ac2(72(crj, aj) = ^^^^2^ and = A — Ai — A2. If we sum all the contributions at order g"^, 



the total result is 





Figure 4: Triple-exchange: (a) 
planar diagram; (b) non-planar di- 
agram. 



W(Ci,C2) 



A2 



2iVM4 



AiA2{AiAs + A2A3 + 3A1A2). (3.7) 



A remark on the Sg2_g2 contribution is in order. This is the 
only contribution in a SU (N) theory and one can verify that its 
small r— expansion is in agreement with the OPE result [ p.ll |, 
supporting the idea that the leading contribution to the Wilson- 
loop is determined only by Oj^. 




Figure 5 

diagram. 



We now come to considering the g^ contribution. Since, 
at this order, the = 4 interactions will start contributing, 
a complete analytic evaluation of all the relevant integrals is 
out of reach. However one can write compact formulas which 



(a) X-diagram; (b) H- 



can be used as a starting point for a numerical evaluation |18]. 
We shall exploit this possibility in a future paper. Here we 
shall instead be interested in singling out the coefficients of 
contributions of the form r'^log(r), potentially present in the 
evaluation of the connected correlator. The knowledge of these 
coefficients already provides non trivial information on the properties of the correlator. In fact, as 
explained in the previous section, a non- vanishing result for these coefficients would clash with the 
expectation that the correlator localizes. 

^Only at this order do the interactions start contributing to the connected Greens functions. 
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For this computation, we limit our attention to the gauge group SU {N) and we can separate the 
diagrams into two classes: the ladder diagrams and the interaction diagrams. The ladder diagrams 
are depicted in fig. ^ and it is easy to realize that they cannot generate any contribution of the 
form r'^log(r). They are actually analytic in the small r— limit. The contributions r'^log(r) are 
instead generated by the interactions diagrams in figs. |5| and ^. The origin of this non analytic 
behavior can be traced back to the small distance singularities appearing in the integration over 
the position of the vertices. Thus in order to extract these logarithmic singularities, we have to 
first perform these integrations analytically, and only after that can we expand in powers of the 
radius. To illustrate the procedure let us start by considering the X-diagram. Its expression can 
be cast into the following compact form 

X = ^ dTidT2daida2 [(xi o y2){x2 o iji)- 

-(xi oX2)(yi oy2)]I^'^\xi,X2,yi,y2), 
where {x o y) = x ■ ij — \x\\y\@x ■ &y with |x|0^ = MjeirsX^x'^ and 



T(^) (xi,X2,yi, 2/2) = 7773^ / 7Z ZXtJZ -TT^. ZJ^TT. IaI"' ^^'^^ 



d'^w 

{2tt)^ J (xi — w)'^{x2 — wY{yi — w)'^{y2 — wY^ 

Here and in the following Xj = x(rj) and yt = y{(Ti) will denote points on the upper and lower 
latitudes respectively (see fig. |l|). The integration over w in ( p.S| ) can be performed and it is then 
straightforward to extract the singular part when we shrink the latitude 9 = 61 lo the north-pole 
of the sphere 5^ (see appendix A for details.) The singular part is given by 

jW-"^-(xi,X2,yi,y2) = -3^x 

^ da ^ 

(1 - a){yi - 0:2)^(^2 - xif - a{l - a){xi ~ X2Y{yi - 1/2)^ + a(yi - xiY{y2 - X2Y ' 



I 



where r = sin^i. The integration over the circuit is straightforward and can be evaluated by 
Taylor-expanding in r. At leading order we find that 

5r^ cos^ ( ^ ) log(r) 



X-"- = ^^li + O(r^). (3.11) 

Consider now the //—diagram in fig. (^). We can write the contribution from this diagram as 
follows 



Ai Bi 



+ Q^\xi,yi,w)D^P^^X2,y2,w) + Q^^xi,yi,w)n,,Q^''{x2,y2,w) 
B2 A2 



(3.12) 



where 

f27r 



P^{xi,yi, w) = I dndai [2yi^ {xi ■ dy^Ii{xi,yi, w)) - 2xf{yi ■ dx,Ti{xi, yi, w))] (3.13) 

and 







f-2n 

Q^' {xi,yi,w) = / dTidai{xioyi){d^Mli{xi,yi,w) - d Mli{xi,yi,w)). (3.14) 
Jo ' 



_ » _ 



In eqs. ( gl3|) and the index M is a ten-dimensional label running; from 1 to 10 and 

in particular we have defined = [x^ ,iQ^\x\) and Bm = (5^,0). The function Ii{xi,yi,w) is 
defined by the scalar integral 



Ii{xi,yi,w) 



(27r)6y {^Xi-zY{y,-zY{w-zy 



(3.15) 



The spatial components of P satisfy the following two simple identities: z^P^ = d^P^ = 0, as 
can easily be checked by direct computation. Moreover, for two latitudes parallel to the plane (2, 3), 
P^ and trivially vanish. Since P'^ is a just a function of z^, all these properties are consistent 
if and only if = 0. This result simplifies dramatically the computation for the correlator of two 
latitudes: in fact the contributions Bi and B2 in ( 3.121) are identically zero. Recall, in fact, that 
Q^^ is different from zero (by construction) only when M is spatial. Thus we are just left with Ai 
and A2 to be computed. 

Let us first compute first A2. It is convenient to rewrite this contribution as follows 

A2 = J dTidT2daida2Xi o yiX2 o y2{dxi - dy^) ■ {dx2 - dy2)n{xi,yi;x2,y2), (3.16) 



where 



W(xi,yi;x2,y2) 



d zd w 



{2Try^ J (xi — zy(yi — z)'^{z — wy{x2 — wY{y2 — 



(3.17) 



The action of {dxi — dy^) ■ {dx2 — dy^) on 'H{xi,yi;x2,y2) can then be evaluated with the identity 
(A. 7) given in [^]. One finds 

{dxi -dy^)-{dx2 - dy2)n{xi,yi;x2,y2) = 

^j(^)(xi,yi,X2,y2)((xi - X2f{yi - y2f - {xi - y2f{x2 - yif')+ 



{xi - yiY{x2 - y2Y L' 
1 



+ 



(y(xi,X2,?/2) - y{yi,X2,y2) + Y{x2,xi,yi) - Y{y2,xi,yi)) 



(3.18) 



where Y{xi,X2,xs) = Ii{xi, X2, xs)[{xi — xs)^ — {xi — X2)'^]. When the first latitude {9 = 9i) is 
shrunk to zero the logarithmically divergent terms can be generated by X*-^^ and by the Y that 
depends both on xi and X2- Therefore we can write 



\3 

» sing. _ ^ 

^2 -8iV2 



, , , {xi o yi){x2 o y2) 
dTidT2daid(T2- 







'{xi - yiY{x2 - y2f 



j(4)--^-(xi,yi,X2,y2)((xi -X2)2(yi -2/2)' 



{xi - y2f{x2 - yif) + 



1 



(2vr)^ 



: {Y^'^s- (a;^ , X2 , y2) + Y^'^a- (^^ , xi , yi)) 



(3.19) 



where we have defined 

^^'"^•(xi,X2,y2) =xr^-(xi,X2,y2)[(xi -y2)' - {xi - X2f] 

and 

Y''^^-{x2,xuy2)^ir'-{xi,X2,y2)[{x2-y2? - {X2-Xif]. 
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The expression for x^*"^' is given in appendix A. The integration over the circuits can then be easily 
performed with the help of Mathematica if we first expand the integrand of ( |3.19| ) in powers of r. 
At leading order we find 

Tr^cos^ (^) log(r) 

^5'°'= — M — +°<'')' 

To complete the evaluation of the //—diagram we have to compute the contribution Ai. The first 
step is to add two total derivatives to the integrand of 



P^'{xuyi,w) = dn da 
Jo Jo 



2yi (ii • dy^Ii{yi - w,xi - w) - Xi ■ dxil2{yi -w,xi- w))- 



2xf{yi ■ dx^Iiixi - w,yi - w) - yi ■ dy^l2{xi - - w))) 



K2 



(3.21) 



These two new terms obviously yield a vanishing result when the integration runs along the circuits. 
The function l2ix, y) is defined in appendix A. Since the following identity for Xi and X2 holds |10] 



^ :Ii{x,y) - 7^X2(2;, y) 



1 X- log ^ 



^x^^^'^'^'^^ ^y^'^^^'^'"' 32tt^ [{x - y)'^ - y^Y 
the combination Ki appearing in can be rearranged in the following compact form 



(3.22) 



Ki=- 



+ 



1 



d 



647r^(yi — wy dri 
1 (xi — w) ■ Xi 



{xi — w)"^ 

{x2 - yif 



+ 



log 



?>2'k'^ {xi — wY {yi — wY \{x2 — J 



[xi - yi) 



(3.23) 



The combination K2 can be also recast into the same form. The only difference from (3.23) is that 
the roles of xi and yi, and of ri and cji, are exchanged. The terms in Ki and K2 that are total 
derivatives with respect to ri and ai can be dropped since they yield a vanishing contribution to 
P^^, and we are left with the compact expression 



P^\x,,yi,w) 



1 



2n 



yi^ (xi - w) ■ xi - xf^ {yi - w) ■ yi ( {xi - yif 

dndai 72 log 7 72 

[xi - wy{yi - w)'' \{x2-y2Y 



Then, if we take into account that 

f.27r 



-Uy,P^'^ {x2,y2,w) = J dTidoi [22/2*^x2 • dy^ - 2xfy2 ■ d^^] 
we can rewrite the Ai contribution in the following form 



(27r)4 (x2 - wf{y2 - wf 



(3.24) 



(3.25) 



^1 = 7779 / dTidT2daida2 log 

Jq \{x2-y2Y 



[{yi o y2)x2 ■ - {yi o X2)y2 ■ 9^2] X 



X xi ■ S{xi,X2,yi,y2) - [(x'l o y2)x2 • dy^ - {xi o X2)y2 • dx2]yi ■ S{xi,X2,yi,y2) 



where 



S'^(xi,X2,yi,y2) = 



d^w 



(xi - wY{yi - wY{x2 - wY{y2 - wY ' 



(3.26) 



(3.27) 
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The nice feature of ( ^.26 ) is the disappearance of one of the integrations over the position of 
the vertices. Although this result simplifies the procedure for extracting the logarithmic terms 
appearing in the limit 6i — > 0, the computation is still a little bit cumbersome and some of the 
details are given in appendix A. Here we shall only give the final result after the integration over 
the circuits. At the leading order in r(= sin^i), we find 



r cos 



log(r) 



5127r2 



(3.28) 



The final set of diagrams to compute are depicted in fig. ^. We have two contributions that we 
call respectively lYup [(c) in fig. |6| and lYdown [(d) in fig. and a diagram which takes into 
account the one-loop correction to the effective propagator [(e) in fig. |6|. We shall denote this third 
diagram by Budiag. To begin with we focus our attention on lYup, whose expression is 



lY 



up 



8A^2 



/ dTidT2dTsda2e{Ti,T2,T3){{xi o y2)x2 ■ {dy2 - dx^) - {±1 o X2)y2 ■ dx^]li{xi,X2,y2), 

Jo 

(3.29) 

and on lYdowm which is obtained from lYup by ex- 
changing the roles of a and r (and therefore Xj and 
yi). Here J is the constant defined by the integral^ 

f27r 

J= I dai{xioyi)D{xi - yi), (3.30) 






where D{x) is the usual Feyman propagator. When 
we shrink the upper circle to a point, the logarithmic 
Figure 6: The two "I Y- diagrams" and the behavior can originate only from lYup. The contribu- 
self-energy correction. ^^^^ lYdown yields analogous behavior when we shrink 

the lower circle. However, when evaluating lYup, we 
also encounter divergences at coincident points (ri T2) in the integration over the upper circuit. 
This singularity though is compensated by the standard ultraviolet-divergence of the self-energy 
graph: half of diagram Budiag cancels the divergence for ri — > T2, while the other half cancels 
the same singularity in lYdown ^or ai — > o"2. Therefore, in order to safely extract the logarithmic 
behavior when we shrink the circuit to zero, we have to first realize this cancellation. 
To begin with, performing a trivial integration by parts, we can rewrite lY^p in the following form 



lY 



up 



/ dTidT2dT3da2e{Ti,T2,Ts) {{xi o 2/2)2x2 • - {xi o ±2)1/2 ■ dx2}Il{xi,X2,y2)+ 
Jo 



2/ dTidT2da2{xioy2)Ii{xi,X2,y2) +- dTidT3da2{xi o y2)Ii{ 

Jo ^ Jo 



x\,xi,y2] 



(3.31) 



The singular part for coincident points is now singled out in the last term, which is proportional 
to Xi(3;i, xi, 7/2)- Since 



Budiag 



8iV2 



dTidTsda2{xi o y2)Ii{xi,xi,y2). 



(3.32) 



This integral is independent of n, namely it is constant, because the integrand is function only of ai — n and we 
are integrating a periodic function over the interval [0, 2tt]. 
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half of Budiag exactly cancels the singularity present in Yup and we are left with 
lY, 



up 



"sFP / dTidT2dTzda2e{Ti,T2,n) {(xi o 7/2)2x2 • dy^ - {xi o ±2)1/2 • dx2}'Ii{xi,X2,y2)- 

-2 dTidT2da2{xioy2)Ii{xi,X2,y2) 

Jo 

(3.33) 

This expression does not exhibit any singularity at coincident points. The logarithmic part arising 
when we shrink the upper circle to a point is then obtained by replacing Ii in the above expression 
with the x^*"^' found in appendix A. Next we Taylor-expand in r and integrate over the circuits. 
At leading order in r we find 

cos^ ( ^ ) log(r) 



Let us now sum all the different contributions at leading order in r 

r^cos^ (f ) log(r) / 5 ^ ^ j \ 

XSing.^jYsmg.^^ sing. sing. ^ \^ 1^ ^ \ g , /ggg) 

7r2 ^768 256 512 15367 

Namely, we have verified that the logarithmic singularities cancel at the first non trivial order. This 
implies that the effective anomalous dimension of the operator Oj^ defined in the previous section 
vanishes at one-loop, supporting the idea that this operator is actually protected. 

As we will show in the next section, this result is consistent with the result coming from the 
zero instanton expansion of QCD2. 

4. The conjectured matrix model description 

In the previous sections we have tried to argue that the correlator of two (or more) Wilson-loops of 
type (|1 . 1| ) might be an exactly solvable quantity since it belongs to a topological sector of = 4. 
This idea, in fact, passes a certain number of non trivial tests: [a] the observable is 1/8 BPS in- 
dependently of the position and the form of the loops P; [b] there is a candidate topological twist 
of the Af = 4: theory, where one of the supercharges preserving the correlator becomes a scalar 
[c] finally, if we compute the behavior of the correlator when one of the circuits shrinks to a point 
we get a smooth limit with no logarithmic singularity. This last property in particular, should be 
contrasted with what happens for the correlator of two circular Maldacena- Wilson loops |l^ : there 
the logarithmically singular behavior was present and signaled the impossibility of a matrix model 



description for this observable |1£]. 



In this section we shall accept this idea, and focus our attention on 
the problem of writing a general formula for the correlator of two 
Wilson-loops. The starting point is to recall that the expectation 
value of one Wilson-loop appears to be computed by the matrix model 
Figure 7: Cylinder ampli- describing the zero-instanton sector of a Wilson loop for QCD2 on the 
tude two sphere 0, ^] . Since the single Wilson loop and the correlator 

generically share the same symmetries we expect that this equivalence 
also extends to the case of correlators. Therefore we conjecture that the correlator of two Wilson 
loops of type (|l.l| ) is given by the multi-matrix model, which evaluates the zero-instanton sector 
of the correlator of two loops for QCD2 on S"^. 
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The construction of this matrix model is quite simple since QCD2 is an almost topological theory 
(it is invariant under area-preserving diffeomorphisms) and its observables can be computed with 
the help of some simple string-like Feynman-rules p6[| . For the present computation we need just 
three ingredients: the cylinder amplitude (heat-kernel propagator), the disc and the Feynman rule 
for the observable, i.e. the Wilson loop. The first quantity is represented in fig. |^ and is given by 



}C{A;Ui,U2) = {U2\e-^\U,) = Y,XR{Ui)xUU2)e- 



-C2{R) 



(4.1) 



R 



where A is the area of the cylinder and the sum runs over all the representations R of U{N). 
The amplitude also depends on the two holonomies Ui and U2 defined on the two borders of the 
cylinder. There is in fact a dual representation for the cylinder amplitude where the sum over 
representations is replaced with a sum over the instanton charges 



N 



exp 



. Y.{4>^-9p^^ + 27Tie,) 



(4.2) 



where {e*^'} and {e*"^'} are the eigenvalues of the matrices Ui and U2 respectively and 



J{9^ 



Jj2sin 

i<j 



The disc is obtained from (|4.lD by choosing one of the two holonomies to be trivial - namely 
equal to the identity. Finally, the insertion of a Wilson loop with winding number n is realized 
by introducing the factor Tr(C/") at the border of the cylinder. The amplitude for the correlator 
of two non-intersecting loops with winding numbers ni and 71-2 is schematically represented in fig. 
1^, and the corresponding expression is given by the following two-matrix integral over the unitary 
matrices: 



K(Ai;Ui,1 




Ui 



Tr(U:) 



K(A3;Ui,U2) 




Tr(U?) 



Figure 8: The string- like Feynman-diagram for the correlator of two Wilson- loops. 



W{Ai,A2) =^ JvUiVU2TiiU^^)TT{U^-')JC{Ar,l,Ui)}C{Ar,Ui,U2)JC{A2-,U2,l) 

P&Sn l,m,s&'E^'' \r,s=l 



(4.3) 



i=l 

2 



{9 A2) 2 -.sfAT-DV.m,- A , n N / 1 



l)(A^-i)E,-, A(^. + 27rm,) exp ( -1^;^^ ^{(^i + 2^m. 

" " i=l 



A being the Vandermonde determinant. The amplitude W(Ai, A2) is related to the true correlator 
by the relation W(^i,^2) = 2^y^{Ai, A2), where Z is the partition function of QCD2 on the 
sphere. We can extend the region of integration over the entire IR^^ by means of the sum over £ 
and m and we can rewrite the above expression as 



N 



/ 1 ^ 1 ^ 1 ^ \ 

X exp (^-^ E - ^ E('^^ - + 2^^^^)' - ^ E -^n • 

(4.4) 

The result (|4.4| ) is the exact amplitude and it contains all instantonic corrections. To single 
out the zero-instanton sector of this amplitude it is sufficient to consider the case where all in- 
stanton numbers vanish. If we introduce the diagonal matrices = diag(0i, . . . , ^at) and 
$ = diag(i?i>i, ...,(/) at), using the Itzykson-Zuber integration formula and defining the hermitian 
matrices Vi = U'^OU and V2 = V<^V-^, we can recast the original integral as the following 
hermitian two matrix model for the correlator of two Wilson loops'' 

WiA,A2) =^ / Z?F,Z)y2e-^^^^^^'^-^^^^^^'^^^^^^^^^^^Tr(e-^^)TV(e-^^) = 

=77^ / Z?liW2e-^^^^^'^-^^^^^^)-^^^«^^-^^)'^TV(e-^^)TY(e-^^), 
CnN^ J 



(4.5) 



where the normalization is chosen to be 



= I i^y.i^y^e-^"^^^^^)-^^^^^^^^^^^^^^^^^ (4.6) 



'^The generalization of this result to the case of n loops is trivial 

where A\, An are the areas enclosed respectively by the first and and last loop (by "enclosed" we mean the region of 
not containing other loops) and Ai^i+\ is the area between the th and (i + 1)— th loop. 
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Actually, in the sector Sj = of (4^), the angular integration can be performed by means of 
an expansion in terms of Hermite polynomials and by exploiting the relation between integrals 
over Hermite polynomials and Laguerre polynomials. Then one finds the following finite A'^ closed 
expression for the connected correlator 

W{Ai,A2)-W{Ai)W{A2) = 

_1 (AiA2(ni+n2)2+A3(nf Ai+n|A2))g2 ^ (^^^^ ^ ^^(^^ ^ ^^^^ _^ ^ ^^^^^ \ 



X 



(Ai(A2+A3)n2+A2(Ai+A3)n2)32 
2A 

^ (fmn^A^A^y-^' (ii - 1)! (g^nlA2{A^ + Ar)\ ^,^„,^ ^5X^1(^3 + ^2) 

«1,«2 = 1 



-j^e 2A X (4.7) 



where A = Ai + A2 + A^^ is the total area of the sphere. For small g this expression can be expanded 
in a power series and one finds 

W{M,A2) - W{A,)W[A2) = 

AiA2{AiA2{nl + nl + mras) + ^3(^1^^ + A2n^))g^nin2 

2^2 



g f Ai^A2iA2 + A3f{2N^ + N)nf A^^ A2HA2 + A3) {2N^ + N) n2nl 
' "^"^ [ tIaW^ + 12l^iV^ 

Ai^A2^ (3^3(^2 + ^3)^2 + Ai (3A3A2 + A2 (4^2 + 1))) ns^nf 



-+ 



UA^N 

Ai^A2HAi + A3){2N^ + N)n2^ni ^1^(^1 + ^3)2 {2N^ + N) 
^ 12^3^2 + 24A3A2 j+^Wj- 

(4.8) 



This result, after decompactifying the sphere, agrees with the perturbative results we have obtained 
up to 0{g^) from Feynman graph calculations using the Mandelstam-Leibbrandt prescription for 
the vector propagator in light-cone coordinates [18|. Let us compare the perturbative result ( [4. 81 ) 
with the actual computation in = 4 done in section 3. After performing the standard redefinition 
g'^ I— > —g'^/A and setting ni = n2 = 1, we find complete agreement up to order g'^. Notice, moreover, 
that the agreement with QCD2 demands the absence of logarithmic singularities when the area of 
one of the loops is small, to all orders in perturbation theory. Our g^ result of sec. 3 is consistent 
with this prediction. 

In order to analyse the large N limit, we can write a simple compact representation for the connected 
correlator in AA = 4 SYM by exploiting a contour representation of the Laguerre polynomials 

\{Ai w2n'^-\- A2A2TL2W1 j 

rr./ . . N TTr, A Nrr./. N "'l"'2 f dwi f dW2 6 A^n,.,^^ J^ J^ 

WiA,A2) - W{M)W{A2) = J^ l j- l,^ ^ -2 , (4.9) 



^^2^21^1 - AiniW2^ 



where Ai = A — Ai and A2 = A — A2. This expression can be computed as an infinite series of 
Bessel functions. We limit our attention to the case ni = n2 = 1 and are actually interested in the 
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normalized correlator, which is given by 



A 



W1W2 iV2^2 



A1A2 k 



lA,A 



'=+1 h 2 



1^2 



\A2A2 



k=l 



A1A2 



(4.10) 



In the next section we will be interested in comparing this result with the strong coupling prediction 
of super-gravity. For this reason, we have to expand the above result for large A. This can easily 
be done by recalling that 



h{z) 
h iz) 



1 + 



(4.11) 



Then the correlator in the strong coupling regime becomes 



^ A A1A2 

W1W2 ^ A^ 



'A1A2 



A1A2 

A1A2 V V A1A2 



+ 



(4.12) 



The first term in the expansion corresponds to the U{\) factor present in U{N) and we shall drop 
it since it is not generally considered in the super-gravity analysis. The first non trivial term which 
can be compared with super-gravity is the second one. 



5. Correlator at strong coupling 



We can also use the AdS/CFT correspondence |28] to compute the correlator of the latitudes at 
strong coupling, in the limit where they are well separated compared to their radii, i.e. in the 
limit that they migrate to opposite poles of the sphere. In this limit the correlator is dominated 
by the exchange of light SUGRA modes between the two worldsheets describing the Wilson loops 



at strong coupling [21, |2^, p^ ]. 

Sometimes, as has been the case for certain chiral primary operators, two point functions with 



the Wilson loop can be computed exactly [29, 27 1 in the gauge theory and succesfully compared at 
strong coupling to a SUGRA calculation of the same quantity. Indeed, by taking the "square-root" 
of the contribution to the correlator of two Wilson loops from a specific SUGRA mode, the two- 
point function of the Wilson loop with the operator dual to that mode is recovered In this 
section we will present a striking agreement between the exchange of certain such SUGRA modes 
and the strong-coupling limit of the QCD2 result ( [4.10 ). In order to prove that the QCD2 result 
truly captures the correlator at strong coupling, cancellations between further SUGRA modes will 



have to be demonstrated. We leave this to a further publication [|18[ 



5.1 An intriguing connection 

There appears to be a rather intimate connection between the QCD2 result presented in section Q 
and the two-point functions of latitude Wilson loops with chiral primary operators built upon the 
scalar field ^3. In the work [27[ it was shown that 



{WOjjx)) 
(W) 



1 

2N 



RsinO 



JXsinO- 



sin( 



sm I 



(5.1) 



where is a latitude Wilson loop at polar angle and 



Ojix) 



TV ($3 + i$4) 



(5.2) 



where x ^ Rsm9 measures the perpendicular distance between the operator and the loop. This 



demonstrates that the matrix model which yields (1.2) also captures two-point functions with those 
CPO's sharing a minimum amount of SUSY with the latitude Wilson loop. 

Let us look then at the contribution of the Oj to the correlator of two latitudes, at polar 
angles ^-nd ^i, taken near opposite poles of the sphere to enforce x S> Rs'mO. Note that 
X = R cos Oq — R cos 9i , we then have 



Oj 



A sin ^0 sin 6i 

J=2 



X 



sin ^0 sin Oi \^ 
(cos Oq — cosOiY J 

Ij(y\ sin^o)/j(VA sin 01 
h (VAsin^o) h (\/Asin0i' 



(5.3) 



This expression is valid strictly at leading order in the large separation limit. The reason for this is 
that (U) ignores quantum corrections between the propagators joining the operator to the Wilson 
loop; this is only valid in the strict large separation limit as shown in ||2^, 27 1. The expression 
( p.3| ) bears a striking resemblance to the QCD2 result ( 4.1C| ). In fact, the only difference lies in the 
factor in round parentheses which is risen to the power J. However, taking the large-separation 
limit of this factor, that difference disappears and (5^) is exactly equal to ( [4.10I ). Thus the QCD2 
result gives, in the large-separation limit, exactly the contribution of the exchange of (|5.2| ). This 



agreement is valid at any value of the coupling, and indeed, in |27] it was shown that at strong 
coupling the result is recovered from supergravity. 

At leading order in weak coupling, this agreement is puzzling for the following reason. It is 
not exactly the operator ( |5.2D which is present in the latitudes' OPE, since there is no coupling to 
$4. Indeed, the calculation of the correlator given in ( 2.11| ) shows that all the operators present in 
the latitude's OPE (2.9) participate in the correlator at this order in A. It is therefore a curious 
coincidence that (|5.2|) produces the same contribution at weak coupling (i.e. J = 2) as the true 
composite operator (|2.9| ) present in the actual OPE. Before addressing this issue further, we present 
a remarkable strong coupling calculation. 

It is interesting to go beyond the strict large-separation limit, and test the QCD2 result ( 4.101 ) 
to higher orders in the shrinking radii of the two latitudes. It turns out that at strong coupling, the 
associated SUGRA calculation giving this information is tractable. In keeping with the intriguing 
connection between the contribution of ( ^.2D to the correlator and the QCD2 result, we begin by 
computing the exchange of the SUGRA modes dual to ( |5.2| ) in an expansion about small latitude 
radii 6*0 and 9i (where the polar angle of the latitude at the south pole is given by tt — 6*1). 

The supergravity modes dual to (|5.2| ) are fluctuations of the RR 5-form as well as the spacetime 



metric. They are by now very well known, and details can be found in |21||22]|32]|29||30|. The 
fluctuations of the metric are 



5gai3 



6 J 



{x)Y\^l), 



2Jg^ps'{x)Yj{n), 



(5.4) 
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where /x, u are AdS^ and a, f5 are S"^ indices. The symbol x indicates coordinates on AdS^ and O 
coordinates on the . The bulk-to-bulk scalar propagator for the field {x) is* 



^W\Fi{J, J -3/2, 2J-3; -41^) 



(5.5) 



where in an AdS^ given by ds^ = (dxg -|- dx?)/xQ, = xo2;o/((2;o — xqY + (^^j — XiY). The full 



details of the calculation will be presented in |18], however it is essentially that found in There, 
the strict large-separation limit was employed by setting the hypergeometric function to 1. Here 
we keep higher terms in the expansion. The results are as follows 



J 



A 



{W{x)W{x)) 
{W{x)){W{x)) ~ 8iV2 



5 • 3 • 26 



22 



9l9l + 9le\ ^ 



26 7 • 33 • 26 



3-27 
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5 • 3 • 23 



+ 0{9 



14 \ 



J 



A 



{W{x)W{x)) 
{W{x)){W{x)) ~ 32iV2 



39l9t ^ elef + el9f ^ H9l9l + 9l9t 



25 5-27 

23 {91 91^ + 91^ 9f) 3^9l9f + 9l9l) 



3-26 
J = 4 : 



72-52 ' 5-25 

{W{x)W{x)) A 



7 • 5 • 33 • 27 



7 • 5 • 25 



+ 0{9 



16\ 



{W{x)){W{x)) 256A^2 



q5 /)5 , Pq + Pq t^i ^0 ^1 + ^0 ^1 , ^0 ^1 



3 • 2 



32-22 



32-24 



The QCD2 result ( [4.10D in the large A hmit is 



{W{x)W{x)) 



{W{x)){W{x)) 



A sin 00 sin^i 
4iV2 



J tan 



tan-'—, 



(5.6) 



(5.7) 



QCD2 J=i 
Ignoring J = 1, we may expand in 9 order- by-order in J: 



J 



J 



{Wix)Wix)) 



{W{x)){W{x)) 
{W{x)W{x)) 



A 



{W{x)){W{x)) 



QCD2 



QCD2 



8iV2 
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32iV2 
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23 (e^0|O + C^i 
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7-5-33-27 
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+ O(0^6) 



'See Pl 



for the definitions of ao and Bj. 
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There is a remarkable matching of highly non-trivial terms between these two calculations! The 
difference between the two calculations sets-in quite late 



A 



(SUGRA-QCD2)^=2 = ^ 
(SUGRA-QCD2)^=3 = ^ 



+ 



+ 



3-2'^ 52-3 5 • 3 • 23 



27 



+ 



29 



+ 



72-52 



7 • 5 • 25 



(SUGRA-QCD2),=4 = ^^ 



24 



(5.9) 



Although we have considered values of J up to J = 4, we expect a similar pattern for arbitrary J. 
5.2 Other modes 

The remarkable agreement displayed in the previous section does not prove that the QCD2 result 
captures the correlator of the latitudes at strong coupling. Beyond the issue of the discrepancy at 
order ^g^f, the catch is that the SUGRA spectrum contains two other modes which couple to the 
string worldsheets and also produce ^q^^ terms, thereby potentially spoiling the agreement with the 
QCD2 result. These are the Kaluza-Klein modes of the NS-NS B-field of type-IIB supergravity, and 



have been described in |32], c.f. their equation (2.48) and what follows it. There is a fluctuation of 
the B-field with both legs in the which is described by a scalar of mass-squared —3 (corresponding 
to a gauge theory operator of protected dimension 3) given by 



a" 



(5.10) 



with k = 1. There is also the fluctuation of the B-field with both legs in the AdS^ portion of the 



(5.11) 



geometry 5Buu, which has been discussed in |31[]. It has the Kaluza-Klein expansion 



5B. 



a%{x)Y\n), 



and the leading k = 1 harmonic corresponds to the following protected dimension 3 operator (where 
A, B are SU(4) indices) 

2i$^^F+ (5.12) 

These contributions must cancel out if the QCD2 result is to hold. Beyond these modes, there are 
also fluctuations of the dilaton, massless vector, and massless tensor which provide contributions 
which lead as 6q9\ and must therefore also find a way to cancel each other, should the QCD2 result 
truly describe the correlator at strong coupling. Indeed this is the reflection at strong coupling 
of the curiosity of the fact that the operators of classical dimension 3 and 4 contributing to the 
correlator at weak coupling seem to have the same effect as replacement by (^]^) (with J = 2). 
The full calculation of these SUGRA modes, and the question of whether or not they cancel, will 
be explored in a companion publication ||T8[| . 

The matrix model result (|1.2| ) contains a rescaled coupling constant A' = Asin2 0. The two 
point function of the latitude with the CPO (^) leads as A' but ends-up as ^/y at strong coupling. 
This explains why in the OPE the operator Tr is weighted by 6'^ but ends-up contributing as 9^ 
at strong coupling. The first descendent of this operator appearing in the OPE of the latitude is 
Tr <I>393<I>3 and comes with weight 6^, thus one would expect its contribution at strong coupling to 
be 9^, potentially explaining why the discrepancy between the QCD2 result and the contribution 
from CPO's built on ^3 sets-in at order 9q9i. 
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A. Appendix 



The integral Xi defined in (3.15), for example, was computed in |]T^] and a useful representation for 
the final result is 



Ti(xi,X2,X3) 



da 



1 



647rVo {y — 



log 



a[(x - yY - y^] + y 
a(l — a)x'^ 



(A.l) 



where x = xi — X2 and y = x^ — X2- The only logarithmic behavior in this integral arises when xi 
and X2 approach the same point xq (namely |x| — > 0), and is given by 



1 



da- 



1 



647rVo (y — axy 
Next we consider the integral 

jW(xi 



log X 



1 



da 



log(xi - X2) 



647rt/o {{x3 - X2) - a{xi - X2)y 



d'^2 



(47r2)4 J [xi - zY{x2 - zY{yi - z)^{y2 - ' 
It is well-known that this integral can be computed in terms of X\ [^] . In fact if we define 



X 



^ _ (xi - y2Y _ {X2 - y2Y ^p. _ (2/1 - V2Y 



{x\ - y2) 



2 ' "^2 



we find 

J*^^)(xi,X2,X3,X4) 

Then 



™2 ™2 ^2 



{X2 - y2Y ' 



X Q 



(yi - y2f 



(A.2) 



(A.3) 



(A.4) 



-2-2-2 

2^X2X3 



T^^>'''3-{xi,X2,yi,y2) 
-I 



{Att'^Y J {xi - zY{x2 - zY{x3 - zY 47r2 

log(xi - X2)^ 



Jl(xi - X2,X3 - X2). (A.5) 



2567r6 



da 



(1 - a){yi ~ X2Y{y2 - xiY - a(l - a)(xi - X2Y{yi - 2/2) + a(yi - xiY{y2 ~ X2) 



2 • 



(A.6) 



For our goals, the most convenient way to compute the integral defined in (|3.27D is to use 
the technique of |^5|, which allows us to reduce the tensor integrals to scalar integrals in higher 
space-time dimensions. We shall perform this reduction in 2u) dimensions and for arbitrary powers 
of the denominators. The final result is very nice and compact 



n 

1=1 



W 



2\a4 



x'^&{uj + l;ai + 5ij) 

(A.7) 



_ on _ 



where 

r-(2t^)/' ^ — IT ^("«) f <f'^w 

(W; ai) = J ^^^^ _ ^)2)ai ((^2 - ((X3 - U;)2)a3 ((2:4 - u;)2)a4 " ^ 

In computing Ai we also need the derivative with respect to of the above expression. After 
some manipulation this derivative can be arranged as follows 

^ r{ai) d f vji'cf^w 



n 



L 47r«i+i dx^ J {{xi - ((x2 - w;)2)«2 ((^3 - i/;)2)«3 ((^4 - u;)2)a4 

4 4 (A-9) 
= S'"'&{u; + l;ai + fe) + 2vr ^ ^ - X2)''6(a; + 2; + % + + fe). 

fe=i j=i 

Finally, the only other ingredient necessary for our calculation is the behavior of the integral 
S{2u;ai) when xi and X2 approach the same point xq. 



c^o ^ _ r - ai) r (a; - as) r(a3)r(a4)r (ai + 02-0;) ((xi - a;2)2)('^-«i-«2) 
2567r'=i r (2w - ai - 02) 



1 + 2 { as— + 04— • ( (a;2 - xq) + {xi - X2) ] + 0{{xi - X2) ) 



(A.IO) 
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